BioData Mining
○ Springer Science and Business Media LLC
Preprints posted in the last 90 days, ranked by how well they match BioData Mining's content profile, based on 22 papers previously published here. The average preprint has a 0.02% match score for this journal, so anything above that is already an above-average fit.
Motta, J. A.; Motta, M. d. M.; Fernandez, C.
Show abstract
In this work, we present a machine learning model for identifying pathogenic DNA variants. The model was learned from the analysis of normal and pathogenic sequences extracted from the ClinVar database (supported by NCBI). This analysis was based on a conceptual semantic model of DNA sequences converted to peptide sequences (amino acid sequences) governed by a well-defined grammar, which allowed us to apply NLP techniques, specifically Part of Speech tagging (POS tagging). Our predictive model was built by combining two techniques: CRF (from the Markov model family), which performs the sequencing, and BiLSTM (a deep learning model) which captures the past and future content of the sequences. The training space was created with the sequences of 105 genes associated with approximately 27,000 pathogenic variants. The model was evaluated using the metrics precision, P-R and ROC curves, AUC, and confusion matrices. Its performance was also compared against five known methods for predicting pathogenic variants. The results show exceptional performance that exceeds expectations and places this new method at the state of the art for predicting pathogenic DNA sequences.
Souza, R. d. O.; Rodrigues, W. F.; Couto, B.; Dos Santos, M. A.
Show abstract
Logistic regression remains a widely used classification method due to its interpretability and computational efficiency, but its direct application to high-dimensional biomedical data is limited when the number of features greatly exceeds the number of samples. In this paper, we propose a reformulated logistic regression framework designed for feature selection and classification in complex high-dimensional settings. The method is evaluated on three biomedical datasets, including scenarios with tens of thousands of attributes and substantially fewer samples. Across these datasets, the proposed approach achieved clear separation between control and disease groups while selecting a compact set of features. Several selected features were consistent with previously reported disease-associated markers, supporting the biological plausibility of the model, while additional selected features suggest potential novel candidates for further investigation. These results indicate that the proposed framework may provide an interpretable and computationally efficient alternative for feature selection in high-dimensional computational biology applications.
Yılmaz, A.; Szydlik, S.; Taheri, G.
Show abstract
BackgroundAdverse drug-drug interactions (DDIs) cause preventable hospitalizations, but exhaustive experimental screening of all drug pairs is infeasible. Many computational predictors rely on SMILES or other molecular representations, limiting their direct applicability to biologics and other non-small-molecule therapeutics. We present a structure-agnostic framework that combines semantic representations derived from Medical Subject Headings (MeSH) with graph-derived topology from a Drug-Target-Protein knowledge graph constructed from DrugBank and UniProt. We further investigate how variation in MeSH annotation depth affects predictive performance. ResultsDrugs are grouped according to their deepest MeSH annotation level (Low, Mid, or Deep), and performance is evaluated across the resulting interaction categories in transductive and inductive settings. The Intermediate ontology scope (Low+Mid) provides the most stable performance, while adding Deep-level terms offers limited and inconsistent benefit. Lightweight topological descriptors are integrated with MeSH features through instance-wise, dimension-specific latent-space gating, using curated reliable-negative pairs for supervision. Fusion improves mean performance over the MeSH-only baseline across all six categories in the transductive setting. Under induction, the clearest gains occur for Low-Low interactions ({Delta}AUROC = 0.056;{Delta} F1 = 0.137) and Low-Mid interactions ({Delta}AUROC = 0.077;{Delta} F1 = 0.114). ConclusionsMeSH annotation depth is associated with systematic variation in DDI prediction performance that aggregate evaluation can obscure. Graph-derived topology is particularly beneficial when ontology annotations are shallow. The framework provides a common, structure-agnostic representation compatible with both small-molecule and biologic therapeutics and supports first-pass DDI prioritization for subsequent expert assessment.
Tuerhanbayi, B.; Fan, X.; Wang, J.; Wan, S.
Show abstract
Acute myocardial infarction (AMI) is one of the leading cardiovascular diseases worldwide and remains a major cause of mortality. Early risk prediction can help clinicians identify high risk patients shortly after admission and support timely monitoring and individualized treatment. Previous AMI risk assessment approaches predominantly rely on a single model structure or fixed feature representation, which may limit their ability to capture diverse risk related patterns and reduce predictive performance. To address these challenges, we propose MetaAMI, a random projection based meta-learning framework for AMI outcome prediction. Specifically, patients were selected based on ICD-9 and ICD-10 diagnostic codes for AMI from the Medical Information Mart for Intensive Care IV (MIMIC-IV) v3.1 database. Features were transformed through multiple random projections, with each random projection generating a distinct lower dimensional feature representation. Subsequently, baseline classifiers were trained on each lower dimension representation to predict in-hospital mortality among AMI patients. The predictions were aggregated to construct an integrated feature representation, which was used as input to a meta-learner architecture. By effectively integrating complementary information from diverse baseline models, the meta-learner refined the decision boundary and enhanced overall predictive performance. Survival analysis and SHapley Additive exPlanations (SHAP) analysis were further performed to evaluate clinical utility and interpret the model predictions. Benchmarking results based on MIMIC-IV dataset suggested that our MetaAMI consistently outperformed all the baseline classifiers across seven evaluation metrics including Accuracy, Area Under the Curve (AUC), F1 Score, G-Measure, Jaccard Index, Youden J, and Matthews Correlation Coefficient (MCC). In addition, feature importance analysis showed clinically relevant predictors of in-hospital mortality. In summary, MetaAMI provides an effective and robust solution for machine learning based AMI risk prediction. We anticipate that the application of MetaAMI will have a positive impact on clinical risk stratification and personalized treatment strategies for AMI.
Shenoy, A.; Zekarias, A.; Viklund, A.; Mitchell, J.; Barrett, J.; Sandberg, L.; Meldau, E.-L.; Taavola-Gustafsson, H.
Show abstract
Background Large Language Models (LLMs) are increasingly explored for pharmacovigilance tasks, including information extraction, case documentation, and single-case causality assessment. However, their ability to support causality assessment at the case series level -- a complex, time-intensive task requiring clinical reasoning across multiple reports -- remains unexplored. Objective To investigate how a large-scale general-purpose LLM can support pharmacovigilance professionals in assessing causality in a case series, and to explore how prompt design influences the quality of the model's reasoning. Methods GPT-4o was used to assess causality for five drug - adverse event combinations, using an adaptation of the Bradford Hill viewpoints for case series assessment. The combinations represented varying drugs and vaccines, adverse events, and case series sizes (5-402 reports). One combination served as a negative control. Structured prompts were iteratively developed and refined using one combination, then applied to all combinations. LLM-generated assessments for each viewpoint were qualitatively evaluated by human annotators for accuracy (precision), and the LLM's coverage of key aspects from the original signal text was assessed for one combination (recall). Results Across all five combinations, annotators agreed with 79-92% of the LLM's output sentences. Full disagreement was consistently low (3-7%), with errors typically involving misinterpretation of complex report details rather than outright fabrication. Prompt design substantially influenced output quality; providing Bradford Hill viewpoint descriptions, including case series data, and adding explicit anti-hallucination instructions improved specificity and grounding. For the recall assessment, 15 of 23 key segments from the original signal text were reflected in the LLM output. The overall summary assessments demonstrated balanced reasoning, correctly distinguishing between positive safety signals and the negative control, and provided a coherent synthesis suitable as a starting point for human assessors. Conclusions LLMs have the potential to generate contextually nuanced and largely accurate preliminary causality assessments of case series aligned with the Bradford Hill viewpoints, with a low but non-zero hallucination rate. These findings support LLMs as a tool to augment, not replace, expert judgment in signal assessment. Future work should address larger and more diverse signal sets, improved evaluation frameworks for generative output, and the integration of pre-computed summary statistics to reduce errors.
Prol-Castelo, G.; Tejada-Lapuerta, A.; Urda-Garcia, B.; Nunez-Carpintero, I.; Garcia-Verellen, E.; Montagud, A.; Valencia, A.; Cirillo, D.
Show abstract
Medulloblastoma is a childhood brain tumor traditionally classified into four molecular subgroups. Recent evidence suggests that Groups 3 and 4 represent a biological continuum rather than distinct entities, a paradigm shift with significant implications for understanding disease biology and treatment strategies. Nevertheless, assessing this hypothesis is challenging mainly due to data scarcity. In this study, we analyze the largest available transcriptomics dataset to provide compelling evidence for the existence of an intermediate subgroup between Groups 3 and 4, characterized by distinct molecular features. To overcome limitations posed by data scarcity, we employ synthetic data generation using a Variational Autoencoder and apply explainability techniques to identify key relationships between gene expression and disease subgroups. Furthermore, by incorporating Machine Learning Fairness approaches, we demonstrate that overlooking this intermediate subgroup can result in treatment disparities. Our findings are further supported by both existing and newly proposed studies using diverse datasets and methodologies, including graphbased analyses and multi-scale simulations, underscoring the robustness and reproducibility of our results. This study demonstrates the potential of synthetic data generation to refine rare disease subtyping and advance our understanding of the underlying biological mechanisms. Keywords: Medulloblastoma, pediatric cancer, representation learning, autoencoder, synthetic data
Xuan, H.; Pasupuleti, R.; Liu, B.; Sun, H.; Zhang, J.; Yao, Z.; Zhong, C.
Show abstract
Bioinformatics software and databases are essential components of modern life science research, yet their mentions in the scientific literature are often inconsistent and difficult to systematically identify at scale. The lack of a comprehensive and up-to-date catalog of bioinformatics resources hinders efforts toward automated biomedical knowledge extraction and streamlined data analysis. Here we present SNAIL, a hybrid named entity recognition framework designed to automatically identify bioinformatics software and database (SW/DB) names from biomedical texts. SNAIL integrates complementary lexical and semantic modeling strategies. The lexical component captures orthographic patterns and contextual cues characteristic of SW/DB names, while the semantic component leverages contextual embeddings generated by transformer-based language models such as SciBERT, combined with an explicit token-masking strategy to enhance entity-focused representations. A large training corpus was constructed automatically through a hybrid pipeline that integrates citation-hinted extraction with large language model-assisted distillation. Evaluation on two independent benchmark datasets and real-world research articles demonstrates that SNAIL substantially outperforms existing approaches, including domain-specific methods such as bioNerDS2 and general-purpose large language models such as ChatGPT, Gemini, Grok and Claude. Applying SNAIL to large-scale literature analysis further reveals distinct journal-level preferences across bioinformatics subfields. These results demonstrate that SNAIL provides an accurate and scalable solution for identifying bioinformatics resources in scientific texts and enables systematic meta-analysis of tool usage and research trends.
Byun, J.; Saha, D.; Han, Y.; Shaw, V. R.; Siminovitch, K.; Amos, C. I.
Show abstract
BackgroundGenome-wide association studies (GWAS) often fail to identify higher-order epistatic interactions that contribute to complex inheritance patterns of traits and diseases. While machine learning (ML) can capture non-linear relationships, extracting interpretable insights from these models remains a challenge. We propose a novel tree-based feature engineering framework that uses Classification and Regression Trees (CART) to explicitly encode high-order interaction decision paths as dummy variables. We investigate three path-based encoding strategies: (i) all decision paths, (ii) leaf-node paths only, and (iii) internal-node paths only. This approach aims to transform complex decision boundaries into discrete features that capture nonlinear interactions that are not readily captured by traditional association models. ResultsThe framework was evaluated using genetic data for ANCA-associated vasculitis (AAV). To manage the high dimensionality of the engineered feature space, we applied a comprehensive suite of ML methods across three tasks: (1) Ensemble Learning (Random Forest, XGBoost, and Gradient Boosting Machine); (2) Decision Tree Analysis (CART); and (3) Regression and Classification Tasks (Regularized Linear Regression/LASSO, Support Vector Machine, and Logistic Regression). Stepwise feature selection and regularization were employed to isolate the most informative interaction patterns. Results indicate that incorporating CART-derived interaction paths--particularly those from high-impact regions of the tree--significantly improves classification accuracy and model interpretability compared to using the original feature space alone. ConclusionsThe proposed framework provides a robust, scalable methodology for identifying high-order genetic interactions. By bridging the gap between the predictive power of ensemble ML and the necessity for mechanistic insight, this approach offers a clearer mapping of the combinatorial genetic processes underlying complex diseases. While applied here to AAV, the method is highly adaptable for exploring the genetic architecture of diverse populations and complex traits.
Honore, A.; Rech, T.; Scrivens, A.; Binotto, I.; Zandvoort, C. S.; van der Staaij, H.; Peck, M.; Zivanovic, S.; Stanworth, S. J.; Hartley, C.; Dame, C.; Deschmann, E.; the Neonatal Transfusion Network,
Show abstract
Background and Objectives: Preterm infants are commonly transfused, yet direct cardiorespiratory effects of red blood cell (RBC) transfusions remain poorly understood. We explored the feasibility of using multicentre electronic health data (EHD) to study such cardiorespiratory responses. Methods: Highly granular routine EHD were collected from preterm infants born <32 weeks gestational age at three European centres. Heart rate, oxygen saturation, and respiratory rate were evaluated 12 hours before and after the RBC transfusion. Results: A total of 321 transfusions in 164 infants were analysed. Overall, there was no significant change in the rate of bradycardia and apnoea following transfusion. Cardiorespiratory parameters varied substantially between infants; e.g. 20% of transfusions were associated with an unexpected, significant increase in heart rate. Respiratory rate and oxygen saturation exhibited similarly heterogenous patterns following transfusion. In sub-group analysis, the proportion of transfusions with increased heart rate was significantly higher within the first two weeks than later (32% vs 13%, p=0.0019). Conclusions: Multicentre EHD extraction allows to identify otherwise masked short-term effects of RBC transfusions on cardiorespiratory parameters, possibly indicating cardiac or pulmonary overload. Such effects may vary with adaptation to anaemia. Analysing EHD may ultimately enable personalized transfusion practice.
Fateh, K.; Yerukala Sathipati, S.
Show abstract
Ovarian cancer is among the deadliest gynecologic malignancies, and its molecular heterogeneity limits accurate prognostic stratification. Although multi-omics approaches have improved predictive modeling, many prioritize predictive performance over biological interpretability, limiting their clinical translation. We developed an interpretable three-stage machine learning framework integrating mRNA, microRNA, DNA methylation, copy number variation, and protein expression data from The Cancer Genome Atlas. Hierarchical feature selection was combined with a weighted ensemble of ElasticNet, ridge regression, support vector regression, XGBoost, and random forest models to estimate overall survival time in patients with ovarian cancer. Multi-omics integration outperformed every single-modality model, achieving a Pearson correlation of 0.752, a concordance index of 0.779, and a mean absolute error of 8.57 months between estimated and observed survival time, compared with 0.48 for the best single modality. The framework identified a 20-biomarker signature dominated by tumor-associated macrophage and complement genes. In an independent survival analysis, VSIG4 and CD163 remained significant after false discovery rate correction, and the signature raised the concordance index over clinical covariates alone from 0.615 to 0.686Enrichment analysis implicated PI3K-Akt, MAPK, focal adhesion, hypoxia, apoptosis, and p53 signaling pathways. This framework couples improved prognostic estimation with biological interpretability supporting multi-omics biomarker discovery in ovarian cancer.
Pucci, F.; Hermans, P.; Tsishyn, M.; Cusato, J.; Rooman, M.
Show abstract
Accurately predicting the effects of pharmacogenomic variants is essential for the development of personalized therapeutic strategies, as genetic variability can influence drug response differently across patients. Here, we assessed several computational approaches using a dataset of pharmacogenomic variants with either clinical annotations or functional characterization by deep mutational scanning, compiled from the literature, with an additional focus on CYP2C9, a clinically relevant drug-metabolizing enzyme. Our results show that, despite recent methodological advances, substantial room for improvement remains. In particular, current methods struggle to distinguish gain-of-function variants associated with increased drug clearance and fast-metabolizer phenotypes from neutral variants, whereas loss-of-function variants that reduce drug clearance are predicted more accurately. The integration of structural and evolutionary information appears to be a key strategy for improving performance, with the coevolution-based StructureDCA method achieving the highest accuracy compared with classical genetic variant-effect predictors and recent deep learning approaches, including the pathogenic-variant predictor AlphaMissense and general protein language model-based methods. Finally, our results indicate that computational models can complement in vitro experiments in clinical variant interpretation, as StructureDCA predictions showed better agreement with clinically annotated phenotypes than large-scale deep mutational scanning data in several cases.
Ravichandrajah, H.; Fischer, A.; Tiago Gomez, A.; Hojeij, R.; Goretzki, S. C.; Felderhoff-Mueser, U.; Park, H.-J.; Kernan, K.; Carcillo, J. A.; Dohna-Schwake, C.; Bruns, N.
Show abstract
Background: Risk adjustment for disease severity in pediatric intensive care research commonly relies on clinical organ dysfunction scores requiring detailed clinical and laboratory information, which is often unavailable in administrative healthcare datasets. We therefore evaluated the feasibility of a coding-based Pediatric Organ Dysfunction Index (PODI) derived from International Classification of Diseases (ICD-10) and Operation and Procedure System (OPS) codes, for approximating sepsis-related organ dysfunction and adjusting for disease severity, using the pediatric Sequential Organ Failure Assessment (pSOFA) score as a reference standard. Methods: In this retrospective single-center cohort study, pediatric sepsis episodes treated between November 2011 and November 2021 were identified. Discrimination for in-hospital mortality and calibration were assessed. Agreement between PODI and pSOFA was quantified using Spearman's rank correlation, and organ-specific agreement using sensitivity, specificity, and predictive values. An expanded PODI incorporating additional ICD-10 and OPS codes was evaluated in sensitivity analyses. Results: A total of 488 pediatric sepsis episodes were included, with an in-hospital mortality of 14.1%. The PODI showed good discrimination for in-hospital mortality (AUC 0.85, 95% CI 0.80-0.89), comparable to the maximum pSOFA (pSOFAmax) (AUC 0.78, 95% CI 0.72-0.83) and superior to pSOFA at sepsis onset (pSOFAonset) (AUC 0.73, 95% CI 0.67-0.80). Agreement between PODI and pSOFA organ-specific components varied considerably across organ systems, with the highest sensitivity to detect pulmonary dysfunction. Correlation between both scores was moderate (0.54 for pSOFAonset and 0.60 for pSOFAmax), indicating that comparable predictive performance does not render the scores interchangeable. The expanded PODI improved organ-level sensitivity for selected components but did not meaningfully improve mortality discrimination. Conclusions: The standard PODI may represent a practical approach to adjust for organ dysfunction and therapy intensity in administrative datasets with ICD-10 coding where clinical and laboratory information is unavailable. Given only moderate agreement with the pSOFA, the PODI should be understood as a covariate for risk adjustment at the group level rather than as a substitute for clinical organ dysfunction scores in individual patients. Further validation and refinement in non-sepsis cohorts are required before broader implementation in large-scale administrative research can be recommended.
Bisaso, K. R.; Kadada, K. R.; Bisaso, K. S.; Ette, E. I.
Show abstract
Background: Parametric time-to-event models require specification of a baseline hazard function, which may influence prediction when the underlying hazard shape is uncertain. This study compared conventional joint longitudinal time-to-event models with mechanistic Multi-Task Logistic Regression, which directly models the survival distribution without selecting a continuous parametric hazard family. Methods: A simulated dataset of 100 individuals with longitudinal sum of longest diameters and event outcomes was analyzed using a shared mechanistic tumor shrinkage regrowth model. Event submodels comprised exponential, Gompertz, Weibull, log-normal, log-logistic, and circadian hazards, mechanistic Multi-Task Logistic Regression, and a hybrid neural-mechanistic extension. All models were estimated jointly using shared patient-specific random effects and longitudinal data. Models were evaluated using longitudinal goodness-of-fit, visual predictive checks, five-fold cross-validated inverse-probability-of-censoring-weighted dynamic area under the curve and Brier scores, integrated Brier score, calibration, and event-interval negative log score. Results: Longitudinal parameter estimates and diagnostics were comparable across models. All conventional hazard models produced identical dynamic area under the curve values within prediction windows, although probabilistic accuracy differed. The log-normal hazard achieved the lowest overall integrated Brier score (0.1928). Mechanistic Multi-Task Logistic Regression achieved the highest later landmark discrimination (area under the curve 0.867 versus 0.798 for all hazard models) and the lowest mean event-interval negative log score (2.362). The hybrid model improved intermediate-landmark discrimination but not overall probabilistic accuracy. Conclusions: Mechanistic Multi-Task Logistic Regression provided competitive joint time-to-event prediction while avoiding baseline hazard-family selection. It represents a practical complementary alternative to parametric hazard modeling, particularly when hazard shape is uncertain and dynamic discrimination is important.
Yelmen, B.; Hofmeister, R. J.; Lutsar, V. K.; Finianos, M.; Stone, B. C.; Joeloo, M.; Krebs, K.; Kivistik, P. A.; Smit, S.; Estonian Biobank Research Team, ; Metspalu, M.; Hudjashov, G.; Milani, L.
Show abstract
Since copy number variations (CNVs) in pharmacogenes can cause significant alterations in drug metabolism, their reliable detection is of high importance both for large-scale studies and personalized medicine. Whole-genome sequencing, and specifically long-read sequencing, is the gold standard for CNV detection. Despite increasing availability of these technologies, genotyping arrays are still widely used as cost-effective alternatives in biobank and clinical settings, yet calling CNVs based on array intensity signals is challenging due to low base pair resolution. In this work, we developed a neural network model, nnCNV, to predict deletions in the CYP2C19 pharmacogene region from array intensity signals. We compared our method to the most widely used algorithm, PennCNV, and demonstrated better performance reaching 100% accuracy in the test dataset. Furthermore, we predicted probe-by-probe CYP2C19 deletion coordinates for all Estonian Biobank samples using nnCNV and PennCNV, and validated these predictions using an identity-by-descent (IBD) sharing method, which also demonstrated superior nnCNV performance. For the deletion samples with conflicting PennCNV and nnCNV predictions, we performed PCR analysis for validation, which showed 97% precision for nnCNV compared to 23% for PennCNV. Finally, we assessed the gradient-based feature importance maps and showed that nnCNV utilizes signal intensity information not only from deletion probes, but also from probes in flanking regions. Our results demonstrate that long-range information, which cannot be utilized by hidden Markov models, can improve CNV calling.
Meyer, D.; Popko, N.; Laub, D.; Schofield, P.; Amariuta, T.; Alexandrov, L. B.; Carter, H.
Show abstract
Genetic feature engineering, used in methods such as transcriptome-wide association study, supports gene-trait association testing by aggregating single variants into gene-level features predictive of expression. To evaluate how different model architectures, LD filtering thresholds, and variant prioritization methods affect expression prediction quality, we trained over 3 million models and evaluated their performance in independent cohorts. Using the best performing models to impute expression and immunotherapy response as an example trait, we found a significant association with the reactive oxygen species pathway (p=0.032). Our model training workflow will support genetic feature engineering towards improved complex trait modeling.
Shokrzadeh, A. J.; Shokrzadeh, P.
Show abstract
Accurate prediction of ex vivo drug sensitivity in acute myeloid leukemia (AML) patients from transcriptomic data is a critical challenge for precision oncology. Existing computational approaches have explored uncertainty quantification in cancer drug response prediction primarily using cell line data, while patient-level AML models typically rely on heuristic confidence measures rather than statistically calibrated uncertainty estimates. Here, we present a framework applying split conformal prediction to patient-level AML drug response modeling using the BeatAML 2.0 cohort. We trained Elastic Net and XGBoost regressors on bulk RNA-seq gene expression profiles from 318 AML patients, analyzing 34,764 patient-drug observations across 122 compounds. Baseline models achieved median Pearson R values of 0.291 (Elastic Net) and 0.281 (XGBoost) across 122 drugs. Wrapping these models with split conformal prediction yielded well-calibrated prediction intervals across three confidence levels: empirical coverages of 81.4%, 90.7%, and 95.5% against nominal targets of 80%, 90%, and 95%, respectively. Analysis of prediction interval widths revealed substantial drug-class-specific uncertainty patterns, with HDAC and BCL-2 inhibitors exhibiting markedly higher uncertainty than MDM2 inhibitors, suggesting a potential association between transcriptomic predictability and drug mechanism of action, although several drug classes were represented by only a small number of compounds. Predictive uncertainty was not significantly associated with ELN2017 molecular risk classification (Kruskal-Wallis p=0.395) or NPM1 mutation status (p=0.788). These results demonstrate that statistically valid uncertainty quantification can be achieved for patient-level AML drug response prediction despite substantial biological heterogeneity. to the best of our knowledge, no published study has applied split conformal prediction to patient-level ex vivo drug sensitivity prediction in the BeatAML cohort, providing a principled alternative to heuristic confidence scoring approaches.
NANDI, S.; Sundararajan, Z.; Subirana-Granes, M.; Espinosa, J. M.; Pividori, M.; Sullivan, K. D.; Galbraith, M. D.; Costello, J.
Show abstract
Down syndrome, caused by trisomy 21, increases the risk of diverse co-occurring conditions. With more than 34,000 related publications indexed in PubMed as of early 2026, keeping pace with this expanding literature is challenging. While general-purpose large language models are widely used for information retrieval, they often rely on broad training data rather than specific evidence. Retrieval-augmented generation (RAG) improves rigor and reliability of responses by linking model outputs to source texts. In research, source texts are peer-reviewed articles. Standard implementations treat all manuscript sections equally, allowing background text to rank as highly as experimental results. To focus model outputs on experimentally supported responses, we developed the T21 Research Assistant, a section-aware RAG system that prioritizes Results sections to ground responses in primary experimental evidence. The system draws exclusively from 1,789 open-access Down syndrome publications from PubMed Central, including 327 NIH INCLUDE-funded studies, and uses a multistage pipeline for query validation, retrieval, reranking, synthesis, and citation verification. Built on NVIDIA Nemotron models, it generates structured, cited responses. Evaluation using expert-curated questions demonstrated strong performance, achieving a BERTScore F1 of 0.712 and recall of 0.758, comparable to or exceeding leading proprietary and open-source models. T21 Research Assistant is available at: https://bioinformatics.cuanschutz.edu/t21-res-assi/
Ravi, A. K.; Gopan, G.; Arumugam, S.; Sethumadhavan, A.; Mani, M.
Show abstract
Abstract Background: The stem cell factor receptor or c-Kit is a type III receptor tyrosine kinase, activated by its ligand Stem cell factor (SCF). Up on activation, c-kit induces signaling pathways that regulates blood cell proliferation, survival, differentiation, and migration. Several studies reported that c-Kit/SCF signaling, contributes to the development and progression of acute myeloid leukemia (AML) in patients. However, the downstream proteins regulated by c-kit activation and their clinical significance in AML remain poorly explored. Methods: Human Acute megakaryoblastic leukemia (Mo7e) cells, were-stimulated with SCF and global protein expression were profiled using two-dimensional gel electrophoresis coupled with MALDI-TOF and LC-MS/MS. Differentially expressed proteins were functionally characterized and validated using patient data from the TCGA-LAML and matched normal data from GTEx, GEO datasets, and quantitative RT-PCR. Their diagnostic and prognostic significance was assessed using ROC, Cox regression, LASSO, Kaplan Meier survival analyses, and a prognostic nomogram model. Results: Proteomic profiling identified 14 differentially expressed proteins in SCF-stimulated Mo7e cells, which are predicted to involved in cytoskeletal organization, protein folding, metabolism, vesicular trafficking, and translational regulation. Transcriptomic analysis of the TCGA-LAML cohort revealed significant dysregulation of CFL1, CCT8, HSP90B1, MDH2, EIF5A, GSN, and TPI1. Integrated ROC, Cox regression, and LASSO analyses identified CFL1, CCT8, and GSN as the most robust prognostic biomarkers associated with poor overall survival in LAML patients. Their expression patterns were validated in independent GEO datasets and by qRT-PCR in SCF stimulated Mo7e cells. Finally, a three-gene nomogram model was developed and validated to predict the overall survival probability of AML patients at 1-, 3-, and 5-year time points. Conclusions: This study identifies CFL1, CCT8, and GSN as key downstream effectors of c-Kit signaling as prognostic biomarkers for AML. These findings provide mechanistic insights into c-Kit-driven leukemogenesis and establish a clinically relevant three-gene signature for AML risk stratification and potential therapeutic targeting.
Choi, L.; McNeer, E.; Beck, C. A.; Neul, J. L.
Show abstract
Bayesian borrowing of external information can improve trial efficiency, particularly in pediatric and rare disease settings where patient populations are limited, but may introduce bias and inflate the Type~I error rate when the trial differs from external studies. Recent U.S. Food and Drug Administration (FDA) draft Bayesian guidance emphasizes careful evaluation of external information, prior specification, and assessment of operating characteristics. This paper compares three meta-analytic-predictive (MAP)-based methods for Bayesian borrowing: the MAP prior, robust MAP (RMAP) prior, and self-adapting mixture (SAM) prior. An adaptive platform trial design in Rett syndrome is used as a case study. Simulation studies evaluate frequentist operating characteristics under varying prior--data conflict, between-study heterogeneity, treatment effects, and clinically significant differences (CSDs) for the SAM prior. The MAP prior achieved the greatest efficiency when external and current data were compatible but exhibited the largest bias under substantial prior--data conflict. The RMAP priors improved robustness through fixed robust-component weights, whereas the SAM prior adaptively adjusted borrowing and was less sensitive to prior--data conflict while retaining efficiency gains when the data were compatible. Although the CSD influenced the degree of adaptive borrowing, as reflected by effective sample size, it had only a modest impact on frequentist operating characteristics. Sensitivity analyses using a skeptical robust component yielded similar qualitative conclusions, while accentuating the differences between the MAP and RMAP priors. These findings provide guidance for evaluating and selecting MAP-based borrowing strategies before trial implementation, particularly in rare disease settings, consistent with current FDA recommendations.
Nehra, N.; Swami, R.; Dadi, D.; Mishra, R.; Sharma, U.; Verma, P.; Sen, M.; Dhruw, N. K.; Jha, A. K.
Show abstract
AO_SCPLOWBSTRACTC_SCPLOWGetting clinical data from different sources to "talk" to each other within the OMOP Common Data Model (CDM) is arguably the most tedious part of multi-center research. While this integration is essential, the transformation process is frequently a manual grind, requiring a rare overlap of deep clinical knowledge and technical expertise. In this paper, we present a framework designed to alleviate some of the burden on the researcher by automating data harmonization through two distinct steps: structural schema mapping and terminological standardization. For the structural piece, we moved away from "black box" logic in favor of a stateful workflow managed by large language models (LLMs) and directed acyclic graphs. By profiling EHR data at the source, our system generates context-aware dictionaries that offer ranked mapping suggestions alongside confidence scores. While our benchmarking showed a 97.5% agreement rate at the schema level and an 84% agreement rate at the value level when compared with human experts, the system appears most effective when treated as a "co-pilot" rather than a total replacement for human oversight. To handle value-level standardization, we implemented a hybrid search strategy that pairs the semantic depth of SapBERT embeddings with the literal precision of fuzzy string matching. By using FAISS for rapid similarity retrieval, the engine attempts to resolve messy or "noisy" clinical descriptions to standard OMOP concepts. This approach seems particularly promising for handling the non-standardized labels that often plague smaller, local datasets. Ultimately, our results suggest that this guided approach can shift the timeline for OHDSI-compliant warehousing from weeks of manual curation to a more manageable and scalable pipeline, potentially lowering the barrier to entry for smaller research teams.